Wave packets and initial conditions in 
quantum cosmology 



S. S. Gousheh and H. R. Sepangi* 

Department of Physics, Shahid Beheshti University, Evin, Tehran 19839, Iran 

February 7, 2008 



Abstract 

We discuss the construction of wave packets resulting from the solutions of a class 
of Wheeler-DeWitt equations in Robertson- Walker type cosmologies. We present an 
ansatz for the initial conditions which leads to a unique determination of the expansion 
coefficients in the construction of the wave packets with probability distributions which, 
in an interesting contrast to some of the earlier works, agree well with all possible classical 
paths. The possible relationship between these initial conditions and signature transition 
in the context of classical cosmology is also discussed. 
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1 Introduction 



The question of the construction and interpretation of wave packets in quantum cosmology has 
been attracting much attention in recent years. Numerous studies have been done to obtain 
a quantum theory for gravity and to understand its connection with classical physics. This 
interest is due to our desire for obtaining the basic dynamical laws for predicting the large scale 
structure of the universe and the determination of the initial conditions from which the universe 
has evolved. The latter is of particular importance since if we accept that the fundamental 
laws of physics are quantum mechanical in nature, then the question of initial conditions is 
embodied in finding the initial condition for the quantum cosmological state of the universe, 
see for example The search for the fundamental dynamical laws has been underway since 
Newtonian times but that for the initial conditions is relatively recent and great efforts have 
been made in this direction. This is, of course, not surprising as the underlying theories for 
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describing the fundamental laws of physics are local quantum field theories above the Plank 
length. However, theories explaining the initial quantum state of the universe are non-local. 
They imply regularities in space that have emerged on large cosmological scales. It is only in 
recent years that the progress in observational cosmology on large enough scales of both space 
and time may present us with the opportunity to look at the viability of the predictions of a 
theory of the initial state of the universe. 

The problem of the relation between quantum cosmology and classical physics is an impor- 
tant one that exists even in simple models . Most authors consider semi-classical approxima- 
tions to the Wheeler-DeWitt (WD) equation and refer to regions in configuration space where 
these solutions are oscillatory or exponentially decaying as representing classically allowed or 
forbidden regions, respectively. These regions are mainly determined by the initial conditions 
for the wave function. Two popular proposals for the initial conditions are the no boundary 
proposal and the spontaneous nucleation from nothing . These proposals have been at- 
tractive to many authors because they lead to some classes of classical solutions represented by 
certain trajectories which posses important features such as predicting an inflationary phase. 

In quantum cosmology, in analogy with ordinary quantum mechanics, one is generally con- 
cerned with the construction of wave functions by the superposition of the 'energy eigenstates' 
which would peak around the classical trajectories in configuration space. However, contrary 
to ordinary quantum mechanics, a parameter describing time is absent in quantum cosmology 
so that the initial conditions would have to be expressed in terms of an intrinsic time param- 
eter, which in the case of the WD equation could be taken as the local scale factor for the 
three geometry [Q. Also, since the sign of the kinetic term for the scale factor is negative, a 
formulation of the Cauchy problem for the WD equation is possible. The existence of such a 
sign is one of the exclusive features of gravity with many other interesting implications. 

The construction of wave packets resulting from the solutions of the WD equation has 
been a feature common to most of the recent research work in quantum cosmology j^, 0, H . In 
particular, in reference the construction of wave packets in a Friedmann universe is presented 
in detail and appropriate boundary conditions are motivated. Generally speaking one of the 
aims of these investigations has been to find wave packets whose probability distributions 
coincide with the classical paths obtained in classical cosmology. In these works, the authors 
usually consider model theories in which a self interacting scalar field is coupled to gravity 
in a Robert son- Walker type universe. The resulting WD equation is often in the form of an 
isotropic oscillator-ghost-oscillator and is separable in the configuration space variables. The 
general solution can thus be written as a sum over the product of simple harmonic oscillator 
wave functions with different frequencies. As usual, the coefficients in the sum are chosen 
according to the initial conditions, which are usually specified by giving the wave function 
at -R = 0, where R is the scale factor. The choice of the initial conditions, or equivalently 
the coefficients in the sum, varies considerably amongst various authors. In |]^, §], the wave 
packets are constructed by summing over the first few terms with the coefficients being chosen 
arbitrarily to be equal. However, in the initial state is chosen to be a Gaussian with the 
appropriate symmetry and the undetermined coefficients being set to zero. The resulting wave 
packets, in general, lack one or more of the following desired properties. Firstly, one expects the 
wave packet describing a physical system to posses a certain degree of smoothness. Secondly, 
there should be a good classical-quantum correspondence, which means that not only the 
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wave packet should be centered around the classical path, but also the crest of the probability 
distribution should coincide as closely as possible with the classical trajectory. Also, to each 
distinct classical path there should correspond a unique wave packet. 

The purpose of this paper is to address and clarify some of these issues. We first argue that 
specifying the wave function at i? = does not furnish a complete set of initial conditions. We 
then suggest an ansatz for the expansion coefficients, or equivalently for the initial conditions, 
which produces wave packets with the desired properties mentioned above. The organization 
of the paper is a follows: In section two we give a short review of some of the models that lead 
to the class of WD equations considered here. Section three deals with the construction of the 
wave packets and in section four we discuss the salient features of this work and conclusions 
are drawn. 



2 Robert son- Walker cosmology with coupled scalar field 

In Robertson- Walker cosmology, one often considers the coupling of a scalar field to gravity. 
The resulting field equations include a 'zero energy constraint'. The WD equation in quantum 
cosmology is the result of the quantization consistent with this constraint. In some simple 
minisuperspace models the solutions to the WD equation describe a system of oscillators with 
unequal or equal frequencies. It would therefore be appropriate at this point to give a brief 
review of some of the minisuperspace models leading to the class of WD equations whose 
solutions are to be discussed here. 

Consider Einstein's field equation coupled to a scalar field 

Rfiu - ^9tiuR = f^T^y{(j)), (1) 

AV-f ^ 0, (2) 

where T^y is the energy-momentum tensor coupling the scalar field to gravity and [/ is a 
scalar potential for the scalar field 0. Hereafter, we shall work in a system of units in which 
K = 1. We parameterize the metric as 

where R{t) is the usual scale factor and k = +1,0, —1 corresponds to a closed, fiat or open 
universe respectively. The field equations resulting from (|I]) and (0) with the metric given by 
(E) can be written as 



Here a dot represents differentiation with respect to time. The choice of is important 
and a particularly interesting choice with three free parameters is ^ 
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f/(0) = A + 



m 



sinh (ad)) + 
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In the above expression A may be identified with a cosmological constant, = d'^U / d(lP'\^=Q 
is a mass squared parameter, 6 is a coupling constant and = |. The Lagrangian giving the 
above equations of motion can be written as 

L = -3Rif + 3kR + i?^ - f/(0)] . (8) 

This Lagrangian can be cast into a simple form for k = describing an oscillator-ghost- 
oscillator system [0, |^. Briefly, this is achieved by using the transformations X = R^/'^ cosh(a0) 
and Y = i?^/^ sinh(a0), which transform the term R^U{(j)) into a quadratic form. Upon using 
a second transformation to eliminate the coupling term in the quadratic form, we arrive at new 
variables u and f , which are linear combinations of X and F, in terms of which the Lagrangian 
takes on the simple form 
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where UJI2 = — 3A/4 -|-m^/2 =1= y/m'^ — 46^/2. The classical solutions for A; 7^ are given in [10 
The corresponding quantum cosmology for = is described by the WD equation written as 

r 92 



Hip{u, v) 



Q2 

+ — + 



uJlv^ /■ il){u, v) = 0. 
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Alternatively, one may obtain a minisuperspace model described by the same WD equation 
but with equal frequencies by considering a closed Robertson- Walker universe with a vanishing 
cosmological constant and containing a conformally coupled scalar field . A different model 
leading to the same WD equation with equal frequencies is obtained by considering a Kaluza- 
Klein cosmology with a negative cosmological constant described by the metric 



g = -df + R\t)- 



+ a\t)dp^ 
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[l + kr^/Af 

where a{t) is the radius of the compactified space |TT[]. The Lagrangian describing the above 
cosmology is 

1 



L = --RaR^ - -R^Rd + -kRa 
2 2 2 

where A is the cosmological constant. By defining uj'^ = 
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'|A and changing the variables as 
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L takes on the form 



L = i 

2 
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It is easily seen that the corresponding quantum cosmology is the same as (|T 
frequencies. We therefore focus attention on solutions to equation (plQl). 
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(14) 
with equal 
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3 Wave packets 



Equation (PUf ) is separable in the minisuperspace variables and a solution can be written as 

1pnun2{u,v) = an^{u)Pn2{v), (15) 

where 

aM = (^Y'^^^^^^e-^-'^ (16) 



Here Hn{x) are Hermite polynomials, and the zero energy condition leads to 

ni + ^) co-i = fr22 + ^) ^0-2, ni,n2 = 0,l,2, (18) 



2/ V 2 

The set {'i/'ni,n2('2^5 '^)} span the zero sector subspace of the Hilbert space of of measurable 
square integrable functions on with an inner product defined in the usual way giving 



tpni,n2{u, V)^n[,n'2{u, v)dudv = 5„i,„;5, 



n2,n2 



That is, the orthonormality and completeness of the basis functions follow from those of the 
Hermite polynomials. 

We can construct a general wave packet as follows 

^{U,V) = J2'An,,n2(^n,{u)f3n2{v), (19) 
ni,n2 

where the prime on the sum indicates summing over all values of rii and n2 satisfying the 
constraint (|18|). As the signs of the kinetic terms in equation (|T0|) indicate, we can take v as 



playing the role of the scale factor and hence the initial condition on %jj is specified by 

ipiu, 0) = Y^"cnian^ (u), (20) 

ni 

where the coefficients c„j are arbitrary and the double prime on the sum indicates that the 
sum runs only over those values of ni satisfying the constraint (piSj) with n2 being even. The 
constraint on the evenness of n2 is due to the fact that /?n2(0) = for odd values of n2- We 
choose the coefficients c„ to be the same as those of the coherent states of a one dimensional 
simple harmonic oscillator, that is 
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where xo is an arbitrary complex number. Comparing the coefficients in equations ( P!U| , pilD one 
finds for even values of 77-2 
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(n2/2)!c„, 



l)-2/2n2!' ^^^^ 



and Arii,ra2 are arbitrary for odd values of n2. This arbitrariness is a consequence of not having 
specified dilj{u,v)/dv\^^Q. The WD equation is a second order PDE, and to have a unique 
solution one has to specify both the wave function and its derivative at a given point. In the 
first reference in the author has chosen A„j^„2 = for odd values of n2, which is equivalent 
to choosing dip{u,v)/dv\^^Q = 0. This might seem to be the only choice when considering 
equation (^). Here we argue that this is not the best or the 'canonical' choice. We have 
chosen Arn^n2 to be nonzero and to have the same functional form for both even and odd 



values of n2. This can be easily done using equation (^). 

The classical paths corresponding to these solutions are the generalized Lissajous ellipsis 
which have the following parametric representation 

u{t) = Mocos(ci;it — 6^0), v(t) = vosm{uj2t), (24) 

where the zero energy condition demands uiUq = U2V0, and 6^0 is an arbitrary phase factor. 
The classical-quantum correspondence is established by the following equation 

Xo = .^uoe''". (25) 

V ^2 

Below we will consider several illustrative examples in order to compare our results with 
some of the previous works. Let us first consider the simplest case which is when ui = 002 = oj. 
Upon inspection of equations (|16D through (|19]), we immediately recognize that, since rii = 
n2 = n and Hn{0) = for odd values of n, firstly An^n are not determined for odd values of n, 
and secondly 



^lJ{u,0) = i^i-u,0), (26) 

(27) 



v=0 



v=0 



Therefore the initial state has to be symmetric and we choose it to be two symmetric Gaussians. 
However this is automatically taken into account by the restrictions imposed on the sum in 
equation ( pO]) and the choice of coefficients given in equation (pT]). Also note that equation 
(^) is redundant due to the properties of the Hermite polynomials. 

Figure 1 Shows the square of the wave packet \ip{u, f )p for |xo| = 12 and Oq = 0, summing 
over both even and odd values of n, and the contour plot of the same wave packet along with 
the classical path superimposed on it. In practice, in order to obtain a reasonable graphical 
representation, the minimum number of terms included in the sum, denoted by n^ax, has to 
be of the order of |xoP- As is apparent from these figures, the crest of the wave packet follows 
the classical path exactly. This is to be compared with figure 1 of the first reference in 
which we reproduce here as figure 2 for the same value of rimax as in figure 1. At this stage, 
due to the symmetry of the circle, the only advantage of our choice of the coefficients seems 
to be that the resulting wave packet is smoother. There is another advantage which will show 
up in the asymmetric case considered below. 

One can easily interpolate between a circle and a line segment simply by changing the value 
of 6q between zero and ±7r/2. Figure 3 shows the wave packet for |xo| = 12 and = 1- This is 
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Figure 1: Left, the square of the wave packet \ip{u,v)\'^ for |xo| = 12 and 9q = 0, cji = ^2 = 
1 with n^nax = 130 and right, the contour plot of the same figure with the classical path 
superimposed as the thick solid line. 




Figure 2: Left, the square of the wave packet \'ip{u,v)\'^ for |xo| = 12 and 9o = 0, ui = uj2 = I 
with rimax = 130, summing only over even values of n and right, the parametric plot of {iplu, v)\'^ 
along its crest showing rapid oscillations. The average frequency of these oscillations increases 
as Umax increases. 
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Figure 3: Left, the square of the wave packet !'?/'(«, f)p for |xo| = 12 and Oq = 1, uji = uj2 = 
1 with rimax = 130 and right, the contour plot of the same figure with the classical path 
superimposed as the thick solid line. 

to be compared with figure 3 of the first reference in where both wave packets with = =tl 
are superimposed due to the choice y4„_„ = for odd values of n. This is not a desired result 
since in order to establish a classical-quantum correspondence we need to have a separate wave 
packet for each distinct classical path. At the classical level, a given set of initial conditions, 
e.g. specifying u(0),f(0), uniquely determines a classical path. At the quantum level the 
initial conditions are specified for example by giving ipiu, 0) and dip{u, I (^'^\v=oi ^^i^ ^i^^ 
uniquely determine a wave packet. Of course the initial conditions have to be consistent with 
the symmetries derived from the general structure of the solutions given by equation (^). In 
the present context, the first condition determines An^n for even values of n, and the second the 
odd coefficients. Therefore the simultaneous appearance of two wave packets corresponding 
to two different classical paths is solely due to choosing dipiu.v) / dv\^^Q = or equivalently 
setting An^n = for odd values of n. This is in contrast to the interpretation given by the 
author in [Q where he relates the simultaneous appearance of two distinct branches to the 
superposition principle, and their only means of separation to the decoherence mechanism. As 
we have shown, for the class of problems discussed here, the use of the ansatz for the initial 
conditions given above would result in the selection of a unique branch without resorting to 
the decoherence mechanism. 

As ^0 =t^/2, the semi-minor axis of the ellipse goes to zero. The results are shown in 
figure 4. This configuration is to be compared with the approximate one obtained by adding 
a few terms with equal coefficients in 0. 

Now we discuss the case with unequal frequency and to be concrete we give an example 
with uJi = 3lj2- Then the constraint equation ( |T^ ) gives n2 = 3ni + 1. Therefore at f = only 
the odd values of rii contribute. Hence we have 



ij{-u, 0) 
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Figure 4: Left, the square of the wave packet !'?/'(«, f)p for |xo| = 6 and 6*0 = tc, ui = 
co'2 = 1 with n^nax = 50 and right, the contour plot of the same figure with the classical path 
superimposed as the thick solid line. 



v=0 



(29) 



v=0 



Again equation (|29| ) is redundant and equation (^8|) is automatically taken into account by 
the restrictions imposed on the sum in equation (pOl) . This together with the choice of the 
coefficients given in equation (BT]) gives two antisymmetric Gaussians. Only the coefficients 



^ni,n2 odd values of rii are determined and those for even values of rii are fixed by the 
prescription given above. Figure 5 shows a plot of f )p for |xo| =6 and 6q = 0. There are 
pronounced oscillations where there is an overlap between different branches of the wave packet 
due to the interference effect. This can be clearly seen from a parametric plot of \ip{u,v)\'^ 
along the classical path similar to figure 2, which we do not show here. In figure 5 these 
oscillations are smoothed out and show up as the bump in the center of the figure. Also note 
that there is a good classical-quantum correspondence. 

Let us next consider the case 6o = tt/2 whose wave packet is shown in figure 6. Note that 
there are pervasive pronounced oscillations due to the overlap. This configuration is to be 
compared with the approximate one obtained by adding a few terms with equal coefficients 
in 0. Moreover, the classical path presented in equation (33) of 0] , describes an open curve 
which seems not to be correct. The correct parametric representation of the classical path 
which is a closed curve is given by equation (|1). 



4 Conclusions 

We have described various classical cosmological models leading to a class of WD equations 
represented by equation (p^Of). The general solution to this equation is a superposition of 
products of oscillator wave functions, as given in equation (|19]). We have argued that there 
exists a set of coefficients which produce wave packets with the following two desired properties. 
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Figure 5: Left, the square of the wave packet for |xo| = 6 and 6q = 0, cji = ?)UJ2 = 

3 with Umax = 50 and right, the contour plot of the same figure with the classical path 
superimposed as the thick solid line. 




Figure 6: Left, the square of the wave packet \ip{u,v)\'^ for |xo| = 6 and 9o = n, ui = 
3uj2 = 3 with Umax = 50 and right, the contour plot of the same figure with the classical path 
superimposed as the thick solid line. 
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First, they posses a certain degree of smoothness, except when there is an overlap between 
different branches of the wave packet, where we expect some interference causing genuine 
pronounced oscillations. Second, there is acceptable classical-quantum correspondence. This 
correspondence consists of three conditions: The wave packet should have compact support 
centered around the classical path, the crest of the wave packet should follow as closely as 
possible the classical path, and to each distinct classical path there should correspond a unique 
wave packet with the above properties. The set of coefficients presented here as an ansatz 
is built from the coefficients of the coherent states as in ordinary quantum mechanics and 
coefficients derived from the functional form of Hn{0) for even values of n. The resulting wave 
packets can be considered as coherent states in quantum cosmology and they posses all the 
aforementioned desired properties, and the classical paths are general Lissajous figures. We 
have shown these features explicitly in several examples. It is worth mentioning that, for the 
class of problems discussed here, the accomplishment of the last condition for classical-quantum 
correspondence by our choice of initial conditions shows that there is no need to employ 
decoherence mechanism to separate out wave packets corresponding to distinct classical paths. 
However, it is generally accepted that decoherence mechanism plays a major role in quantum- 
to-classical transition, see for reviews. Although we have not shown the results explicitly 



here, we can report that as |xo| decreases the degree of classical-quantum correspondence 
diminishes. Also note that there is a difference between these coherent states and those of 
quantum mechanics in that the former are static. 

We have also shown that there exists a continuous parameter — 7r/2 < 6q < 7r/2, whose 
variation gives a continuous interpolation between wave packets or classical paths correspond- 
ing to different initial conditions. When 6'o = we have the most symmetric configurations, 
i.e. loops which are symmetric about the line u = 0. When 6q = ±7r/2 we have the most 
asymmetric configuration, i.e. curve segments. 

This work is also related to the topic of initial conditions in the context of signature 
transition. For example, at the classical level there is a controversy on the value of R and 



at the hypersurface of signature transition. The authors of reference [0 argue that these 
quantities need to be only continuous across the hypersurface, while in |n] it is insisted that 
in addition to the requirement of continuity, the values of both these quantities should be zero 
at the hypersurface. This choice corresponds to the most asymmetric case 6q = ±vr/2, which 
gives curve segments as classical paths. We speculate that these are unlikely configurations 
due to their pronounced pervasive oscillations for ti;i/ti;2 > 1- 
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